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Abstract 

We use QCD sum rules to obtain the weak parity-violating pion-nucleon coupling 
constant f n NN- We find that f n NN ~ 2 x 10~ 8 , about an order of magnitude smaller than 
the "best estimates" based on quark models. This result follows from the cancellation 
between perturbative and nonperturbative QCD processes not found in quark models, 
but explicit in the QCD sum rule method. Our result is consistent with the experimental 
upper limit found from 18 F parity-violating measurements. 
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In this Letter, we use the method of QCD sum rules with the electroweak and QCD 
Lagrangians to predict the weak parity- violating (PV) pion-nucleon coupling constant, 
fnNN- The theoretical prediction of f n NN is an important and challenging problem. To- 
date, the most accurate PV experiments have only shown 1,2 ) that the upper limit for the 
magnitude of this coupling constant is 3-5 times smaller than the "best value" predicted 
by DDH 3 ) on the basis of a quark model and somewhat smaller than that in a similar 
calculation carried out more recently. 4 ^ Since that time others have tried to estimate f n NN 
by means of chiral soliton models 5 ' 6 -* and QCD sum rules. 7 ) This coupling is of particular 
interest because of its sensitivity to the neutral current contribution of weak nonleptonic 
processes at low energies. ^ 

QCD sum rules have been shown to be able to reproduce known properties of the 
nucleon, e.g., fx p , /i n , gA, and of other hadrons. 8 ) However, they have rarely (if ever) been 
used to predict unknown properties. Keeping terms in the operator product expansion 
(OPE) up to dimension 5, we show that there are two main terms in the sum rule for 
fnNN- the unit operator and a dimension D=3 susceptibility. By using an analogous sum 
rule for the strong coupling constant, g^NN, to evaluate this susceptibility, we are able 
to determine the weak coupling f n NN- An important aspect of the present work is that 
we demonstrate that there is a cancellation between perturbative and nonperturbative 
QCD modifications of the weak process. 

We employ a two point function for the nucleon in an external pionic field. Our 
current is the usual one 9 ) 

Vp (x) = e abc K T ( a ;)C7^ 6 (x)]7 5 7^ c (^), 

fj P (y) = e abc [u b (y) lu Cu aT (y)]d c (y)Yl 5 , (1) 

where e abc is the antisymmetric tensor, C is the charge conjugation operator, and a, b, c 
are color indices. The neutron currents are similar, with the interchange of d <-> u. 
Since the most general weak PV n-N coupling is 2,3 ' 10 ) 

H PV (nNN) = ^p(r x 0^ , (2) 
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only charged pions can be emitted or absorbed. For definiteness, we consider the absorp- 
tion of a 7r + so that an initial neutron is converted to a proton, and the correlator we 
consider is 

n = ij d A xe tx - p < 0|T[77 P (x)77 n (0)]|0 >^+ . (3) 

The general form of II for the parity-violating pion-nucleon coupling, as dictated by 
relativistic invariance, is 

u PV = n e 1 + u oP , (4) 

with p = 7^p M . 

The phenomenological evaluation of the correlator is carried out by deriving a disper- 
sion relation for II through the insertion of a complete set of physical intermediate states 
of spin | in the expression of Eq. 3. Using the usual terminology, we refer to this as the 
right-hand side (RHS). We only use the sum rule for Il e , since the sum rule for n o is not 
as stable. One finds for the parity-violating part of Eqs. 3,4: 

II e (p ) = _ + continuum. (5) 

M is the nucleon mass; and the parameter is related to the amplitude for finding 
three quarks in a nucleon at one point and has been determined in a number of sum-rule 
calculations. 8 '' The double pole term, corresponding to the insertion of the one-nucleon 
intermediate state in Eq. (3), has contributions both from the weak pion-nucleon ver- 
tex and the parity violation in the nucleon state itself. As will be shown below, in 
our microscopic calculation using the two-point form only Z -quark loops in the nucleon 
correlator give the parity-violating vertex correction. As is usual in the method, the 
physical property of interest, f n NN, is obtained by treating the double-pole term explic- 
itly, while the continuum and excited states are included in the numerical analysis via a 
parameterization, as discussed below. 

The microscopic evaluation of II, based on QCD and electroweak theory (the so- 
called left-hand sides (LHS) of the QCD sum rules for II), is obtained by means of a 
Wilson coefficient expansion in inverse powers of p 2 . In this work we keep diagrams up 



to dimension D = 5. The lowest dimensional diagrams which we consider are shown 
in Fig. 1. The higher dimensional diagrams which we include are obtained from those 
shown in Fig. 1 by the substitution of Figs. 2b and c for the pion-quark vertex, Fig. 2a. 
The propagators in coordinate space corresponding to the three diagrams of Fig. 2 are: 

• — * — * 

St = -^r-ng 7TqX7T <qq>6 ab 1 5 , 

Sic = m o < m > 9* q T ■ , (6) 

with X-nQ-wq^j < QQ >=< qiTjjsq > n and < qq > itj =< q i^^gcTjcr ■ Gq > w . Here 
g nq is the pion-quark coupling, which is not explicitly used in the present calculation, 
and G represents the gluon field. The susceptibility Xn enters in the evaluation of both 
strong and weak pion-nucleon coupling constants, while enters only for the weak one. 
We will discuss the treatment of these parameters below. We only consider the even 
sum rule, namely that for IT e ; that for n o involves further unknown susceptibilities. The 
evaluation of the diagrams is straightforward. 

For the weak Hamiltonian, we take H w = ^| ( J^jj + jVJVj) with 

= fr^(l - 75 )dcos0 c 
N» = uY(A u + B ul 5 )u + d^(A d + B dl 5 )d, (7) 

where 9c is the Cabibbo angle and A u , A d , B u , B d , are given by 

1 8 
A u = -(1 - - sin 2 ^), 



1 4 
A d = -3(1- 3 sin 2 ^) 



B u = ~B d = — - , (8) 

with 9 W the Weinberg angle. This is the standard model Hamiltonian, which we use for 
the main part of the calculation. We then discuss the QCD effects on our results. 

Since momentum can be transferred in the weak point-like interaction, shown by wavy 
lines representing Z° in the figures, there is an additional integral to be carried out in 



the evaluation of EL For example, we obtain for Fig. la 

Ii\ a = -2 6 G F sin 2 d w g nq J d D hd D k 2 d D k 3 [h ■ (p - k 2 - k 3 ){p -k x - k 3 )k 2 
+ e -{2(k 2 ■ (p - k 3 )h(p - k) + 2(p -h- k 3 ) -(p-k 2 - k 3 )k 2 k 
+ -3(p - h - k 3 )k 2 h(P - h - k 3 )}} 

[{2-nf D klklkl{p -h x - k 2 f{p -k 2 - h) 2 }- 1 , (9) 

where D = 4 - e is the dimension. There is no PV contribution from Figs, (lc) and (Id), 

and the sum of Figs, (lb) and (le) vanish. The integrals in Eq. (9) are evaluated by 

standard Feynman techniques, with dimensional regularization. The result is 

rriat 2\ G F sin 2 Q w g 1 , q 6 2 1 15 3 

n e (P ) = ^27^ P M-P )(- + y ~ 7p) • ( 10 ) 

We regularize the diagram using mass, vertex, and pion-quark vertex counter terms, 
leading to the one-loop corrections to our diagram shown in Fig. 3. The lowest dimension 
pion-quark vertex and mass renormalization diagrams for f n NN are shown in Figs. 3a-c. 
In our approximation of a contact weak interaction, the contribution of Figs. 3a-c vanish 
under a Borel transformation. The mechanism of Figs. 3d and 3e do not appear in 
the external field method. The only nonvanishing diagrams in the infinite Z-mass limit 
are those shown in Figs. 3f and 3g. With a minimal subtraction scheme we obtain an 
additional composite current, which we call r]v'- 



with ki = p — k 2 — k 3 and q = k 2 + k 3 . This current is used for the vertex regularization 
shown in Figs. 3f and 3g. These vertex corrections give the contribution 

u e { v)(p ) = — — p ln (~P )(; + y - 7) • (12) 

Combining Eqs.(ll0|,12) and taking the Borel transform one obtains for the regularized 
diagram la 

n^)(P 2 )= GF -yyJ g " (f -7)Mf. (13) 
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where Mb is the Borel mass. The other diagrams can be evaluated in the same manner. 
The results from the processes of Figs. 1 and those obtained from Figs. 1 with the 
substitution of Figs. 2b and c for Fig. 2a, with the counter terms given by corresponding 
substitution in Fig 3, are 

+ | Xw aL-^M 2 B E 2 + X -ml a E 1 L~ 4 / 9 ] 



M 2 



= WA^e- A ^(2^-l) (14) 

where a = — (27r) 2 < qq > and = (27r) 4 \ 2 v /g nq . We do not include gluon condensate 
diagrams for f n NN] they are of the same order or smaller than uncertainties of our 
calculation. The factors containing L, L = 0.621 ln(lOMe), give the evolution in Q 2 
arising from the anamolous dimensions, and the Ei(M B ) functions take into account 
excited states to ensure the proper large-Mg behavior. The last line in Eq. (|T4"D is the 
Borel transform of the double-pole term from the phenomenological (right-hand) side, 
Eq. (^) . The direct proportionality to sin 2 9w should be noted. 

Finally, by explicit calculation or Fierz reordering, we can show that the contribu- 
tion for W ± exchanges vanish. Thus, as required by symmetries 3 ' 10 '', we find no charged 
current contribution to the weak PV pion-nucleon vertex; such a contribution requires 
strangeness-changing currents and would thus be reduced by sin 2 Q c <m 0.05. Since we 
neglect strangeness in the nucleon and strangeness-changing currents, we obtain no con- 
tribution. 

As we shall demonstrate below, the first two terms in the theoretical form for Il e 
given in Eq. fll4]) are of opposite sign and tend to cancel. This is a crucial point. For 
this reason it is essential to either determine the value of the susceptibility x-k from g n NN 
or to eliminate it from our equations. We do both as an aid in determining the stability 
of our solutions. First, we determine x-n- directly in terms of g^NN [as a function of the 
Borel mass] by using the sum rule for the strong coupling, which is analogous to Eq. 



(Ill), and attempt to use the result to determine f^NN- Second, we eliminate x-k from 



the PV and strong coupling sum rules and find that we can determine f W NN in terms of 
QttNN ■ Details are given below. 

We use the correlator given by the two-point function of Eq. (|3|) for the strong as 
well as the weak interaction. The general form differs from Eq.(^) by the presence of a 
7 5 in each term. The phenomenological (RHS) for the strong pion-nucleon coupling is 
now given by 

U s( p 2 )RH s = ^ + continuum. (15) 

(p 2 — M 2 ) 

Unlike the weak PV pion-nucleon coupling, the evaluation of the strong one leads to 
a problem in that there is no double pole on the right-hand (dispersion relation) side. 
However, as shown by Reinders et. al. u \ the value of the coupling constant g n NN found 
in this way is virtually the same as that found by means of a 3-point function, which 
circumvents the lack of a double pole problem. 

Keeping terms up to D=6, shown in Fig. 4, for the theoretical side (LHS), and taking 
the Borel transform we obtain the sum rule for the strong pion-nucleon coupling: 

9,nn\ 2 n e~ M2/M - = M%L~^E 2 - M% X ^E X + -^L^ 

+ < > E Q M% _ < glQ2 > EoMl{ U _ lnM% (l6) 

where < g 2 c G 2 > is the gluonic condensate. 

Before we discuss our detailed evaluation of the sum rules to obtain our estimate of 
fnNN, let us discuss the structure of Eqs. (|b|) and (|T6|). First, as we discuss below, if we 
use the method of Ref. (12) [which uses arguments of PCAC within the sum rule context] 
to evaluate Xn we find that Xti-Q=-88 GeV 2 . With this value, the Xn term dominates both 
Eqs. ( |1^ , |T6|) with the result that g^NN — 155 [in contrast to the experimental value of 
13.5]. With this value of Xn we find that f n NN > 10~ 6 , at least an order of magnitude 
larger than experiment. 

Secondly, since Xn is the only unknown in Eq. (0), we can estimate the vacuum 
susceptibility using the experimental value of g^NN = 13.5: this gives x-k & — -1.88 GeV 2 , 
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two orders of magnitude smaller than the value given by the method of Ref. (12) [see 
discussion below]. With this value one finds that the first two terms in Eq. Ql4]) , the 
leading terms for f n NN, almost cancel. Note that the second term involving Xn enters 
with the opposite sign in the two equations for f^NN an d g^NN, respectively. This is 
the source of the very small parity-violating pion-nucleon coupling in comparison with 
quark model: there is a cancellation between the dimension zero model-like term using 
perturbative quark propagators and the vacuum pion susceptibility term. 

However, we find that the sum rule obtained for f^NN [Eq.flHD], using the value of 
Xtt(M^) extracted from Eq. (0), is not stable in Mg. Therefore we cannot obtain a 
reliable estimate of f^NN by this method. 

We find that we can obtain a satisfactory sum rule to determine f n NN by eliminating 
from both Eq. (14]) and Eq. (|16|) by taking derivatives with respect to M B . With 
this procedure, and taking the ratio of the weak to the strong sum rule we obtain the 
new sum rule for the weak in terms of the strong coupling constant: 

g.NN Cw N (2Mf, + 3M B -9M 2 N M B )(12E 2 M B + 3<G 2 >E o y 1 ' 

where c w = GpsirfQwQj — 7)/(247r 2 ) = 5.5 x 10~ 8 GeV~ 2 . The sum rule is quite 
stable with a plateau in M B in the region expected, as shown in Fig. 5. Because of the 



strong cancellation between the first two terms in Eq. (FJ) [dimension and dimension 
2 terms], the dimension four term with the unknown parameter mj is important for the 
final numerical value of f n NN- We have taken hiq = in Fig 5. Guided by the value of 
the parameter m needed in the nucleon sum rule 8 -*, we evaluate the sum rule given in 
Eq. 15 with nig taken over the range 0.0 to +0.8 GeV. From this procedure we find: 

U NN « (1.9 to 2.4) x 10- 8 /or 
ml = (0 to 0.8)GeV. (18) 

For negative values of mj the value of f n NN becomes smaller and even negative, but we 
did not find stable solutions for sizable negative values of this unknown parameter. To 
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be consistent with the neglect of gluon condensate terms we quote as our central value 
of fnNN that with tUq = 0, shown in Fig. 5, as 

f nNN « 1.9 x 1(T 8 . (19) 

This coupling constant is an order of magnitude smaller than the "best values" of Refs. 3 
and 4. As emphasized earlier, this result follows from the cancellation of the two leading 
terms in Eq. fll4]) . The first LHS term in that equation, a unit dimension term which 
would correspond to a quark model type calculation, gives a value for f n NN ~ 2 x 1CT 7 , 
similar to the quark model value. The second term, involving the nonperturbative QCD 
vacuum susceptibility, Xn, strongly cancels the first term. Because of this cancellation, 
we cannot expect Eq. fllTf ) to be very accurate, but we find a clear explanation for the 
small value of f w NN, consistent with experiment. 1 ' 2 ) 



The results given in Eqs. [18] and [19] have been obtained using the Hamiltonian of 
the standard model (see Eqs. |7| and |[). Let us now consider the strong interaction 
modifications. These have been estimated in Refs. 3 and 4 using the renormalization 
group method. In the notation of Ref. 3, the operators involved in our calculation are O4 
and O5. Using the tables in Refs. 3 and 4 we find that the our parameter for the parity 
violation, AaB u — A u Bd, would be changed by less than a factor of two in magnitude. 
Since the same parameter appears in all terms, this gives the overall uncertainty arising 
from strong interaction modifications. Therefore, the main conclusion of our work is not 
changed. 

There are two relevant features that we would like to point out. The first one is that 
the use of pseudovector coupling also circumvents the problem of a lack of double pole 
for the strong interaction constant. For the Lagrangian 

CirNN = — $N ll^T ■ TpNVpfa (20) 

we can treat as a constant external axial vector field. The QCD sum rule is then 

identical to our calculation of gA- 8 ^ At the quark level, we have 

Arw = TT$* '/MVA , (21) 



where /„- is the pion decay constant. From our previous result for g A ^ we then obtain 

q'ttNN _ 9a 

, 2M_g A M 
9-kNN — 9-knn —i — \ lz ) 

TY^n Jn 

which is just the Goldberger-Treiman relation. 

As a second feature we wish to attempt an independent estimate of Xn- For this 
purpose we first use PCAC to obtain 

_ f rrj* 

< 0\u i-f 5 u - d i-f 5 d\7T° > = J ^ £ e~ iq - x 

\/2m q 

= x *<t>n < qq > e~ iq ' x , (23) 

where we take the 7r-quark coupling to be unity in this discussion. We then use the work 
of Belyaev and Kogan 12 ) , which assumes saturation of a sum by one pion states: 



< 0\q i"f 5 T 3 q\0 > 7r = J d*xe lQ x < 0\u i^u - d i^dln >< n\q i-y 5 T 3 q\0 > Q ^ 

<qq> (24) 
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2m2 

As described above, the value of x-k obtained in this manner is more than an order 
of magnitude larger than that found by using the value of g n NN from experiment. Once 
more we point out that if we use it in Eq. ( |I6|) we find an order of magnitude discrepancy 
with the strong coupling constant, g n NN- Furthermore, it is clear that this value of Xn 



is inconsistent with Eq. (14), since by eliminating it with derivatives with respect to 
the Borel mass we obtain results an order of magnitude different than with its use. We 
conclude that Eq. (p4|) cannot be correct. We are not certain where the method of 
Belyaev and Kogan errs, but we believe that it is suspect. Note that Xn = (fn/m q )Xn ~ 

20Xn 

In conclusion, we find that the weak PV pion-nucleon coupling due to neutral cur- 
rents is as small as that due to charged currents, ~ 2 x 1CT 8 . This result agrees with the 
conclusion of the chiral soliton model of Kaiser and Meissner 5 ) , but not that of Kaplan 
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and Savage 6 ^. Our result also disagrees with quark model calculations 3 ' 4 ) and with a 
previous QCD sum rule calculation. 7 ) If the coupling is as small as we estimate, it cannot 
be separated from the charged current contribution and thus cannot be found experi- 
mentally; and it is unlikely that the anapole will be seen. 13 ) Although we have omitted 
gluon condensate corrections to the PV correlator, our result is sufficiently small that 
these corrections will not alter our conclusion. Finally, we point out that in the two-point 
QCD sum rule method used here, the small value of f n NN which we obtained is the result 
of a cancellation between a process which can be treated in quark models and a vacuum 
process identified in the method of QCD sum rules. 
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Figure Captions 

Fig. 1. Lowest dimension quark diagrams for the PV weak pion-nucleon vertex. The 
dashed line represents a charged pion and the wavy line a Z°. 

Fig. 2. Quark propagator modifications in an external pion field. 

Fig. 3. Pion-nucleon weak vertex correction diagrams. 

Fig. 4. Diagrams contributing to the calculation of g n NN- 

Fig. 5. Solution for f w NN from Eq. 15 as a function of M B . 
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